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ABSTRACT 

Why do galactic bars rotate with high pattern speeds, when dynamical friction should 
rapidly couple the bar to the massive, slowly rotating dark halo? This longstanding 
paradox may be resolved by considering the dynamical interactions between the galac- 
tic disc and structures in the dark halo. Dynamical friction between small-scale halo 
structure and the disc spins up and flattens the inner halo, thereby quenching the dy- 
namical friction exerted by the halo on the bar; at the same time the halo heats and 
thickens the disc, perhaps forming a rapidly rotating bulge. Two possible candidates 
for the required halo structures are massive black holes and tidal streamers from dis- 
rupted precursor halos. More generally, gravitational scattering from phase-wrapped 
inhomogencitics represents a novel relaxation process in stellar systems, intermediate 
between violent relaxation and two-body relaxation, which can isotropize the distri- 
bution function at radii where two-body relaxation is not effective. 

Key words: galaxies: formation - galaxies: evolution - galaxies: haloes - galaxies: 
kinematics and dynamics. 
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1 INTRODUCTION 

Most of the mass in disc galaxies resides in three compo- 
nents: a thin, rapidly rotating disc of stars and gas, a cen- 
trally concentrated stellar bulge or spheroid, and an invisible 
dark halo that extends far beyond the outer edge of the stel- 
lar disc. The bulge rotates less rapidly than the disc, and the 
rotational properties of the halo are unknown. The halo is 
believed to form by hierarchical clustering of dark matter, 
which results from gravitational instability of small irregu- 
larities present at much earlier times; its structure depends 
on the c osmological model, but usually the halo is stro ngly 
triaxial (Dubinski & Carlberg 1991; Warren et al. 1992) and 
rotates only slowly (median spin parameter A ~ 0.05; see 
for example Steinmetz & Bartelmann (1995). The disc is 
believed to form from baryonic material that collects i n the 
potential well of the dark halo . Numerical simulations ( Katz 
& Gunn 199l|; iDubinski 1994 kepner & Bryan 19981) show 



that the growth of the disc modifies the structure of the in- 
ner halo (by which we mean the part at radii smaller than 
the outer edge of the disc) : the halo becomes more nearly ax- 
isymmetric (b/a <: 0.8) and oblate (c/a ~ 0.6) but remains 
slowly rotating. 



In this paper we examine the long-term gravitational 
interactions between the disc and inner halo. In particular 
we shall discuss the rate and consequences of angular mo- 
mentum transfer from the disc to the halo, which tends to 
flatten and spin up the inner halo, while causing the disc 
material to thicken, heat, and drift towards the centre. Any 
irregularities in either the disc or halo tend to contribute to 
angular momentum transfer; we shall argue that the trans- 
fer is probably sufficiently strong to transform the inner halo 
from a slowly rotating oblate spheroid supported by pressure 
anisotropy into a rapidly rotating thickened disc. 

Rapid rotation of the inner halo may resolve a persis- 
tent puzzle about the structure of barred disc galaxies. A 
rotating bar embedded in a non-rotating dark halo loses an- 
gular momentum to the halo through dynamical friction, on 
a time-scale that is generally much less than a Hubble time. 
Thus bars should have small pattern speeds, corresponding 
to large values of the dimensionless ratio 1Z = (corotation ra- 
dius) / (radial extent of the bar)-which is normally ~ 1 when 
a bar forms. This argument was first explicitly discussed by 
Weinberg (1985), although earlier numerical models of disc 
galaxies with live halos by Sellwood (1980) exhibited rapid 
angular momentum transfer from the disc to the halo af- 
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ter a bar formed in the disc. Hernquist & Weinberg (1992) 
simulated systems containing a non-rotating halo and a rigid 
bar (with no disc) and concluded that the gravitational cou- 
pling was 'sufficiently strong to remove all of the bar's angu- 
lar momentum on a time-scale much shorter than a Hubble 
time'. Little & Carlberg (1991a,b) followed the long-term 
evolution of barred disc galaxies with live halos using a two- 
dimensional code, which greatly improves spatial resolution 
but requires a somewhat unrealistic flat halo. They found 
that TZ increased from 1 to 2 over ~ 30 initial bar rotation 
periods. The most recent and comprehensive simulations of 
this process are by Debattista & Sellwood (1998), who fol- 
lowed self-consistent iV-body models of a barred disc and 
initially spherical halo. In the simulation with the most mas- 
sive halo, they found that the bar pattern speed dropped by 
a factor of five (to TZ ~ 2.3) within 50 initial bar rotation 
periods. Only in models with low-density halos did the bar 
pattern speed remain high. 

Thus a wide range of analytic arguments and numerical 
simulations concur that high-density halos strongly deceler- 
ate bars. On the other hand, observations show that most 
bars are rapidly rotating. Bar pattern speeds can be esti- 
mated in several ways: (i) By comparing the shocks in mod- 
els of gas flow with dust lanes in barred spiral galaxies. Using 
this approach Athanassoula (1992) estimated TZ = 1.2±0.2. 
(ii) From the rotation field of the old disc population and 
the equation of continuity. This method has been appl ied to 
the SB O galaxy NGC 936 and yields TZ ~ 1 .4 ± 0.3 jKent 
& Glau |dell 1989] |Merrifield fc Kuijken 1995] ). (iii) By iden- 
tifying photometric and kinematic properties of stars and 
gas with Lindblad or corotation resonances; this approach 
is probably less reliable since some disc features may be 
transitory and also the pattern speed of the spiral struc- 
ture ma y differ from that of the bar (Sellwood & Sparke 
1988^ 

recent review by Elmegreen (1996) concludes that 
TZ — 1.2 ±0.2 in early-type galaxies; bars in late-type galax- 
ies may rotate more slowly but so far there are very few 
observations. 

Thus, there is a significant and growing contradiction 
between theory (TZ 1) and observations (TZ ~ 1). Debat- 
tista & Sellwood (1998) argue that the absence of slowly 
rotating bars can only be explained if the halo contributes 
less than ~ 25% of the radial force in the region where bars 
are found (around 2 disc scale lengths); this in turn requires 
that the disc provides almost 90% of the total circular speed 
(7 = «disc/wtot > 0.87; see Oiling 1995 for a more precise def- 
inition). The relative contribution of the disc and dark halo 
to the radial force in sp iral galaxies is controversial and un- 



mal dark halo-have a median 7 of only 0.85 (Sackett 1997) 



-7? 

certain ( Bottcma 1997|; Brocils & Courteau 1997 ;_ dourteau 



& Rix }997| ; |5ackctt 1997| ; |Fuchs 199$ |Sellwood 1998| ), but 
it seems unlikely that the very high value of 7 required by 
Debattista and Sellwood is present in most disc galaxies. 
For example, a careful recent discussion of mass models for 
the Galaxy is given by Oiling & Merrifield (1998); the nine 
models they describe with various values for the solar circu- 
lar speed and radius have 7 in the range 0.66-0.83 and their 
favored model has 7 = 0.72. Even 'maximal disc' models of 
external galaxies-which have the largest possible disc mass 
consistent with the observed rotation curve and an isother- 



which is slightly lower than Debattista and Sellwood require. 

A rapidly rotating inner halo could resolve this apparent 
contradiction, by reducing or even reversing the drag on a 
rotating bar*. 



2 DISC-HALO GRAVITATIONAL TORQUES 

Disc-halo angular momentum transfer can result from struc- 
ture in either the disc or halo. We first examine the effects 



of disc structure (§2.1) and find that no known disc features 



transfer angular momentum rapidly enough to spin up the 
inner halo. Halo structures (6E.2I) are more efficient. 



2.1 Structure in the disc 

The discs of spiral galaxies exhibit a wide variety of struc- 
ture, including molecular clouds, spiral arms, and central 
bars. All of these can transfer angular momentum to the 
halo. 

To estimate the transfer rate we adopt a simple model 
in which the circular speed of the disc, v c , is independent 
of radius, and a fraction fa of the gravitational force in the 
disc plane is due to a spherical halo; the remaining fraction 
fd = l — fh arises from the self-gravity of the disc. In this 
model the halo de nsity and disc surface density are (Binney 
fe Tremaine 1987|) 



Ph(r) = f h 



S d (r) = f d 



2%Gr' 



(1) 



If the torque per unit disc area exerted on the halo is N, 
then the disc material in the radius range [r, r + dr] loses 
angular momentum at the rate 2irrNdr; if we assume that 
this angular momentum is gained by the halo material in the 
same radius range, with mass A-kt 2 ph(r)dr and characteristic 
specific angular momentum ~ rv c , then the halo is expected 
to flatten on a characteristic time 

_ 2p h (r)r 2 v c = f h vl 
T ~ N 2-kGN' ( ' 

First consider the effects of discrete masses in the disc. 
We suppose that a fraction /a/ of the disc mass is in objects 
of mass M and radius R. Each such object loses specific 
angular momentum to the halo by dynamical friction, at a 
rate (Binney & Tremaine 1987, eq. 7-24) 



-0.43/ h — In A, 
r 



(3) 



where A « r/R. The torque per unit area on the halo is then 
N = fM^d\L\, and the halo flattening time is 



2.3- 



v c r 



f M f d GM\nA 



* Debattista & Sellwood (1998) find that halo rotation is not 
sufficient to suppress bar slowdown, but their rotating halos are 
constructed by reversing the angular momenta of stars in a spher- 
ical halo, and thus do not rotate as rapidly as a flattened halo. 
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The largest disc objects are giant molecular clouds (GMCs). 
These have a wide range of masses but their contribu- 
tion to dynamical friction depends on the weighted average 
(M 2 ) I (M), which is dominated by the largest clouds at mass 
roughly 10 6 M (e.g. Combes 1991). The largest GMCs may 
not be bound, but unbound objects transfer angular momen- 
tum just as well. The fraction of the total mass in molecular 
gas ranges from 4% for Sa galaxies to 25% for Scd galax- 
ies (Young & Scoville 1991) so fiufd = 0.1 is a reasonable 
average. We then see from equation (^) that GMCs do not 
transfer significant angular momentum to the bulk of the 
inner halo (although individual GMCs at radii ^ 1 kpc do 
spiral inwards in less than a Hubble time; see Stark et al. 
1991). 



and then merge to build up successively more massive but 
less dense objects. The process terminates at a redshift z 
given roughly by 1 + z ~ fi.^ 1 , where 57 m is the density 
parameter. When a small halo merges with a larger one, it 
loses energy by dynamical friction and spirals towards the 
centre until it is tidally stripped or completely disrupted, at 
a radius r ~ r s (M/M a ) 1 ^ 3 , where r a and M s are the radius 
and mass of the small halo and M is the mass of the larger 
halo interior to r. The disrupted halo is stretched into a tidal 
streamer; the length of the streamer after a time t is roughly 
A<£ ~ (t/P) (r s /r) ~ (t/P)(M a /M) 1/3 , w here P is the or- 



bital period (Tremaine 1993; Johnston 1998). Such stream- 
ers can survive as distinct-though unbound-structures in 
phase space for a Hubble time or longer; since the halo ma- 
terial is collisionless the streamers never overlap in phase 
space but simply become longer as phase mixing proceeds, 
forming a kind of phase-space spaghetti. Structure resulting 
from incomplete phase-mixing is already well-known in sev- 
eral other contexts: moving groups in local disc st ars (|Eggen 



Npft wp eaamiBfi aagukt mameataim ijaaafci te soi. 



ral stnjcjmrs in the disc. A tightly wrapped spiral pa.t- 



196!: ; Dehnen 1998), shells in elliptical galaxies (Hernquisi 
Uumn 198$ [Malm 199Sj), clumps in kinematic surv eys of 



tern with azimuthal wavenumber m, pattern speed Q p , ra- 
dial wavenumber k and fractional surface-densit y amplitude 
Ei /E d exerts a torque per unit area on the halo ( Mark 1976 ) 



N = 



->7/2 3/2 „,2 „ ^2^2 

2 7T m ilpphG E d 2 
■ (Ijl/ljd) ■ 



3 fcV 3 v-v-", . (5) 

where a is the one-dimensional velocity dispersion in the 
halo and the halo distribution function is assumed to be 
Maxwellian. In our model a = Vc/2 1 ^ 2 , and using equations 
(jjj) and (Q) the halo flattening time is 

37T 1 / 2 m 2 4 /E d \ 2 

r = —nJ! {cotl) IetJ ; (6) 

where i — tan -1 (m/kr) is the pitch angle of the pattern. 
An optimistic set of parameters is m — 2, i — 20°, / d — 0.7, 
ttp 1 ~ 5xl0 7 yr, andEi/E d ~ 0.2, which yields r ~ 10 12 yr, 
too long to be of interest. Unusually strong and open spiral 
patterns could transfer angular momentum more rapidly, 
but these are normally transitory, induced for example by 
an encounter with a nearby galaxy. 

Angular momentum transfer by dynamical friction on a 
central bar has been discussed in the Introduction. Although 
the frictional forces are strong enough to despin the bar, it 
is much harder to spin up the halo, because the reservoir of 
angular momentum in the bar is too small. Bars typically 
contain ^ 30% of the disc luminosity and extend only to 
~ 30% of the disc radius, so their moment of inertia is only 
a few percent of the disc's; thus transferring all of the bar's 
angular momentum to the halo-as noted, a rapid process- 
would still impart negligible halo rotation. 

We conclude that the known disc structures cannot spin 
up the inner halo. 

2.2 Structure in the halo 

The dark halo is believed to form by gravitational instability 
of small density fluctuations in the early Universe. In stan- 
dard hierarchical models, dense, low-mass halos form first, 



the metal- weak halo ( Majewski, Munn fc Hawlcy 1996), a nd 
the disrupted Sagittarius dwarf galaxy (Ibata et al. 1997). 

A rotating disc exerts dynamical friction on a tidal 
streamer. To estimate the frictional force, consider the ide- 
alized case in which a single streamer of length L & r pierces 
the disc at right angles. At distances <C L, we can approxi- 
mate the streamer as an infinite, straight wire, and we can 
approximate the wire as stationary since the velocity vector 
of its material is approximately aligned with the streamer. 
We can approximate the disk as a uniform, infinite, zero- 
thickness sheet moving past the streamer at velocity v c . If 
the linear density in the streamer is A = M s /L, the frictional 
force exerted on it from disc material with impact parameter 
b < L is 

_ 4^ 2 G 2 A 2 S d b 1 

-Tdrag — , , (I ) 



where &i is the maximum impact parameter considered. This 
drag force is exerted only while the streamer intersects the 
disc. If the streamer follows a circular orbit of radius r, then 
it intersects the disc for a fraction L/(2nr) of its orbit, so 
the average drag force is 

_ 27rG 2 A 2 £ d Lfei _ 27rG 2 M s 2 S d bi 

\-Tdrag) — n — o 



L ' 



(8) 



the maximum impact parameter is b 1 ~ L. Disc material 
with impact parameter 6 > L sees the streamer as a point 
mass, and the time-average of the frictional force component 
from this material is 

2 i/2 G 2 M 2 SdlnA 2 1/2 G 2 A 2 E d L 2 lnA 



(F, 



drag/ — 



(9) 



where A = 62 /&i , and 62 and bi are the maximum and min- 
imum impact parameters considered. Equations (||) and Q 
are similar except for the Coulomb logarithm in the latter, 
but the ratio 62/61 is not usually large, so material with 
b <C L and 6 3> L contributes comparable drag forces. 

Estimating the frictional force is more delicate for long 
streamers (A<f> < 2w). A long streamer intersects the disc 
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several (JV ~ A<f>/2n) times. If the locations of these inter- 
sections are uncorrelated, then the total drag force would be 
given by the drag per streamer (eq. |^ with 61 ~ r) times JV, 



of mass M (e.g. blac k holes). Then the torque on the halo 
per unit disc area is (Ostriker 1983) 



JV = 2 5/2 Ji In A 



G 2 ME dPh r 



(15) 



(■Fdrag) 



2-kG 1 ! M s 2 £ d 
v 2 rA<j> 



(10) 



If on the other hand the intersections are localized in the 
disc-as they would be if the streamer had a planar orbit 
of small or moderate eccentricity-then the JV streamers of 
linear density A act like a single streamer of density JVA. 
Replacing A by JVA and bi by r in equation Iph we find 



where I\ = 0.474. Using equations (^) and (^]), the charac- 
teristic flattening time for the halo is then (cf. eq. GJ) 



iw \ — G 2 M%£d 

(,-Fdrag,) - 2 



(11) 



2 3 / 2 h InA/d GM 

'l0 6 M Q \ / 10 



3kpc 



We regard equation ( |ll| ) as more reliable than ( ]10| ) in most 
cases. 

Equations (Q) or (Q) for short streamers, and equation 
( pd[ ) for long streamers, can be summarized by 



(-fdrag) 



M 



( — ) 



(16) 



Observational upper limits to the mass of compact halo ob- 
jects are reviewed by Carr (1994). Masses M^3x 1O 6 M 
are strongly disf avored because they cause excessive heating 
of the local disc (Lacey & Ostriker 1985; Lacey 1991 ), exces- 



where g ie a dimoneionloee conetant of ordor unityi Tho fric 

tional tbiqiie flattens the 'fapampr nrhit intn trip Hicr planp 

in a time 



gravi tational lensing of gamma-ray bursts (|Marani et al 



sive distortion of gravita ti onally lensed radio jets (Wambs- 
ganss & Paczyhski 1992; |Garrett et al . 1994|), unobserved 



199S ) , a nd excessive growth of the central black hole in the 



Galaxy (Ku fc Ostriker 1994) 



M s v c 



(13) 



fdrag gG 2 T,dM a ' 

For a numerical estimate of the flattening time, let f 3 
be the ratio of the mass of a streamer to the halo mass inside 
radius r, f s = GM S / fhvjr, and use equation ([]]) to eliminate 
Ed. Then equation ( |l3| ) yields 

2n r 



gfsfhfd v, 
3 x 10 10 



yr 



0.3 



3 kpc 



fh 
200 km 



0.01 



(14) 



Streamers containing even a few percent of the the halo mass 
can be dragged into flattened orbits within a Hubble time. 

The principal effect of the frictional torque from the disc 
is to increase the z-component of the angular momentum of 
the streamer. The orbital energy of the streamer changes 
more slowly: since the gravitational field of the streamer is 
approximately stationary, the energy of both the streamer 
and disc are conserved. The total angular momentum of 
the streamer also changes relatively slowly: for example, the 
torque on an axisymmetric streamer is always perpendicu- 
lar to the streamer plane and hence perpendicular to the 
angular momentum vector, so the magnitude of the angu- 
lar momentum is unaffected by the friction. This process of 
increasing L z at fixed L and E rotates the streamer planes 
into the disc plane, without strongly affecting their shape. 
At the same time the disc stars lose L z at constant E, which 
heats t he disc without strongly affecting the radial mass dis- 
tributk n. 



Some authors have argued for much more stringent up- 
per bounds on the mass of compact halo objects, but we find 
these arguments unconvincing, (i) Rix & Lake (1993) esti- 
mate that halo objects with M > 6 x 10 3 Mq would cause 
excessive heating of the stars in the dwarf irregular galaxy 
GR 8. This is a very low- luminosity galaxy [Mb = —10.6, 
corresponding to 2.5 x 10 6 Lq or 2 x 10 _4 L*), so even a dark 
halo containing ~ 30 times the mass in visible stars will com- 
prise only JV ~ 10 halo objects of 3 x 1O 6 M0. The crossing 
time in the halo is t c ~ 10 8 yr, and in Nt c ~ 10 9 yr the halo 
will evaporate, leaving a single or binary halo object at the 
centre, which of course would not heat the stellar distribu- 
tion. This argument predicts that the dark mass in GR8 may 
be concentrated at its centre, producing a quasi-Keplerian 
circular-speed curve, but testing this prediction is difficult: 
in fact GR8 is one of the few galaxies with a declining HI 
rotation curve but th ere is a large and uncertain correctio n 
for pressure support (Carignan, Beaulieu & Freeman 1990). 
(ii) Moore (1993) and Klessen & Burkcrt (1996) argue that 
massive halo objects disrupt globular clusters; the second 
of these papers, which is more thorough and conservative, 
sets an upper limit M < 5 x lO 4 Af0. This result is based 
on the assumption that the disruption time td for a par- 
ticular subset of globular clusters (central density between 
10 3 and 1O 4 M0 pc -3 ) must be long enough that more than 
1% of them survive; for a Poisson process with destruction 
probability per unit time t^ 1 this requires that ta > Q-22tn 
where tu is the age of the cluster population. However, if 
the present clusters represent the survivors of a much larger 
original population-as indicated by the efficiency of the de- 
struction mec hanisms associated with known components 
of the Galaxy ( Gnedin & Ostriker 1997 ; Murali & Weinberg 



To close this section we examine the rate of angular mo- 
mentum transfer if the halo is composed of compact objects 



1997)-then we expect a wide range of survival times in any 



snapshot of the population. In particular, if the distribution 
of lifetimes at formation is a power law, then we expect the 
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fraction of clusters with lifetimes less than td "C tu to be 
~ td/ts so ~ 20% of clusters would be expected to violate 
the Klessen-Burkert criterion. 

We conclude from equation ( |l6| ) that if the halo is com- 
posed of 10 6 M Q black holes, well below the allowed upper 
limit of 3 x 10 6 Mq, then we expect substantial halo flatten- 
ing and rotation out to ~ 1 kpc, which may be sufficient to 



ett & Pogge (1995), who found much flatter halos with axis 
ratios 0.3-0.5. 

High-resolution HI data can be used to determine the 
thickness of the HI layer in highly inclined disc galaxies; 
at and beyond the edge of the optical disc this thickness is 
sensitive to the shape of the dark halo. The most th orough 



suppres s halo dynamical friction on rapidly rotating bars. 1996| ), yields c/a ~ 0.2^'^. Oiling & Merrifield (1997) de 



analy sis of this kind, for the Scd galaxy NGC 4244 (Oiling 

+0.3 

o.i 



3 CONSTRAINTS ON HALO SHAPES 

The arguments in the preceding sections suggest that the in- 
ner parts of dark halos may be thick, rapidly rotating discs. 
As noted earlier, bars are strongly coupled to the inner halo 

and tin's srirmlrl also rntat.p rapidly as nbsprvprl Wp now 



ask wh ether other observations constrain the shapes of in- 
ner halos. Although halos formed by hierarchical cluster- 
ing are generally triaxial, the subsequent growth of the disc 
tends to make the inner halo more nearly axisymmetric (cf. 
§1), and this expectation is confirmed by the upper limits 
to the ellipticities of most dis c galaxies induced by the halo 



Kui- 



potential-typically ^ . 1; see (Pranx fc de Zeeuw 1992 
jken & Tremaine 1994 ; Rix 1996] ). Thus we assume that the 
inner halo is axisymmetric, and restrict our attention to the 
flattening of axisymmetric halos. 

Rix (1996) and Sackett (1996, 1998) review the tools 
that have been used to constrain the halo shape in our own 
and other galaxies, including the shapes of X-ray halos, the 
kinematics of polar-ring galaxies, the thickness of HI layers, 
and gravitational lensing. 

The kinematics and spatial distribution of Population 
II objects can be used to co nstrain the shape and mass of 
the d ark halo in our Galaxy (Monet, Richstone & Schechtei 



1981). Unfortunately the best available results are model- 



dependent and have large error bars: 0.3 < c/a < 0.6 from 
Binney, May & Ostriker (1987) and c/a > 0.34 from van der 
Marel (1991). 

The ellipticity of the mass distribution in edge-on disc 
galaxies can be determined from the shape of their X-ray 
isophotes. This method yields c/a = 0.60 ± 0.13 for the SO 
galaxy NGC 1332 (Buote & Canizares 1996); this is greater 
than expected for disc-like dark matter, but the uncertain- 
ties are large and the result is dominated by halo radii well 
beyond the outer edge of the disc, where the flattening mech- 
anism discussed here would not operate. 

Steiman-Cameron, Kormendy & Durisen (1992) fit a 
precessing-disc model to the dust lanes in the SO galaxy 
NGC 4753 and derive 0.84 < c/a < 0.9; however, this result 
depends on an assumed age for the event that warped the 
galaxy disc. 

The axis ratios of the dark halos in polar-ring galax- 
ies have been estimated by several investigators. Whitmore, 
McElroy & Schweizer (1987) examined three SO galaxies 
with polar rings and found axis ratios for the potentials 
of 0.86 to 1.05 (±0.2), corresponding to axis ratios for the 
density of roughly 0.58 to 1.15 (±0.6); however, two of these 
galaxies were re-examined by Sackett et al. (1994) and Sack- 



rive c/a = 0.75±0.25 for the Galaxy using both the thickness 
of the HI layer and local stellar kinematics. 

Gravitational lensing by disc galaxies can constrain 
the halo mass and shape. Koopmans, de Bruyn & Jack- 
son (1998) find c/a > 0.5 for the edge-on lensing galaxy 
B1600+434. Flattening the halo does not significantly affect 
the overall lensing cross-section of a galaxy but does en- 
hance the fraction of images that exhi bit characteristic ge- 



omet ries associated with disc lensing (Keeton & Kochanek 



199S ) . Unfortunately, disc galaxies are expected to comprise 
only 10-20% of gravitational lenses so a large sample will be 
difficult to obtain. 

Evidently the present data allow a wide range of halo 
shapes, with some indication that moderately flattened halos 
(c/a ~ 0.5) are more common than spherical or disc-shaped 
halos. The data provide no strong evidence for or against ro- 
tating thickened discs of dark matter such as those proposed 
in this paper. 



4 DISCUSSION 

Structure in galactic discs and dark halos couples these two 
components together, leading to angular momentum trans- 
fer from disc to halo, spinup and flattening of the inner halo, 
and heating and thickening of the disc. The time-scale for 
this process is highly uncertain. Known structures in the 
disc are probably not able to flatten and spin up the inner 
halo significantly in a Hubble time, but halo structure is 
likely to be more efficient. In particular, we have identified 
two types of halo structure that may be able to spin up the 
halo: massive black holes, and tidal streamers resulting from 
hierarchical halo formation. 

A rotating flattened halo may be needed to explain the 
high pattern speed of bars (§[j]). Flattening the inner halo 
also thickens and heats the inner disc; this process may form 
some galactic bul ges, which, like discs, are meta l-rich and 
rapidly rotating (Wyse, Gilmore & Franx 1997). Another 



way to make bulges from discs is thro ugh a buckling in- 



stabi li ty of bar structu res in the discs (Combes & Sander: 
1981; Raha et al. 1991). The hypothesis that some bulges 



may be produced from discs is supported by the inner cut- 
offs observ ed in disc surface -brightness profiles as they enter 
the bulge ( Kormendy 1977), and the similar scale lengths o f 
inner discs and bulges (Courteau, de Jong & Broeils 1996). 
Kormendy (1993) has argued persuasively that some bulges 
are really discs in terms of their dynamics and origin. On the 
other hand, the centres of bulges have much higher phase- 
space density than discs and s o cannot arise from discs 
through a collisionless process flCarlberg 1986] ); moreover 
the strong similarity between bulges and elliptical galaxies 
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of similar luminosity suggests that most bulges are formed 
in a manner similar to the family of equivalent ellipticals. 

Dynamical friction on tidal streamers provides a natural 
mechanism to produce flattened dark halos of non-baryonic 
material, and thus offers a counter-example to the usual be- 
lief that if the dark halo is disclike it must be baryonic. If this 
process is important then we expect the inner halos of spiral 
galaxies to be substantially flatter than halos of ellipticals. 

Relaxation from tidal streamers can also be important 
in elliptical galaxies. Gravitational scattering of stars by 
streamers can enhance the relaxation rate compared to the 



usual C handrasekhar formulae for two-body relaxation (Bin- 
ney & r >emaine 1987) and hence isotropize the distribution 



function at radii where two-body relaxation is ineffective 
(energy relaxation is less effective than isotropization be- 
cause the potential from a long streamer only varies slowly 
in time). 

We may contrast relaxation from tidal streamers with 
two other relaxation processes in stellar systems: two-body 
and violent relaxation. The two-body relaxation rate is de- 
termined by assuming that the stellar system is as smooth 
as possible, so the only potential fluctuations arise from 
Poisson noise due to individual stars. Violent relaxation 



(Lynden-Bell 1967) is modelled by assuming that the po- 



tential fluctuations in the stellar system are as large and 
rapid as possible (spatial and time scales of order the sys- 
tem size and dynamical time), an assumption that is only 
valid d uring the initial collapse of the system. Relaxation 



precisely into the disc midplane. However, if there are several 
axisymmetric streamers with the same radius and different 
midplanes, they cannot all be dragged into the disc mid- 
plane since this would lead to an increase in the phase-space 
density of the halo material, which violates the collisionless 
Boltzmann equation. 

There are also direct observational tests, (i) Massive 
black holes scatter a fraction of the stars in the bulge and 
inner disc out to much larger radii; these should appear as a 
distinct population of ~ 10 s M@ in metal-rich high- velocity 
stars on nearly radial orbits (tidal streamers are less effec- 
tive scatterers because their potential is softer). Similarly, 
black holes scatter local disc stars to large epicyclic ener- 
gies, c reating a power- law tail containing ~ 1% of the disc 
mass ( Lacey fc Ostriker 1985| ). (ii) Tidal streamers heat 
disc stars by removing angular momentum at constant en- 
ergy, while massive black holes induce a random walk in 
both energy and angular momentum; these processes may 
leave distinct signatures in the phase-space distribution of 
inner disc and bulge stars, (iii) Heating by tidal streamers 
is most effective early in the galaxy's history, while heating 
by black holes continues throughout its lifetime; this distinc- 
tion may be reflected in the age or metallicity distribution 
of bulge stars, (iv) The presence of a thin disc of old stars 
near the centre of a galaxy would imply that substantial 
angular-momentum transfer from disc to halo has not oc- 
curred since the disc was formed. Some ea rly-type galaxies 



contain discs with scale le n gths & 1 kp 



van den Bosch 



from t jdal streamers is intermediate-stronger than two- Jagg fc van der Marel 1998; ; |Scorza et al. 199% if these are 



body relaxation, and weaker but longer lasting than violent 
relaxation-and arises because a stellar system only phase 
mixes gradually over many dynamical times. Perhaps the 
process should be called 'non-violent relaxation'. 

One important unresolved question is how and when 
the streamers are finally mixed together. Small-scale irreg- 
ularities in the disc and halo (GMCs, globular clusters, etc.) 
can mix the streamers, but relaxation arising from interac- 
tions between streamers may be more effective. This issue 
is also important for experiments that hope to detect the 
phase-space structure of dark matter particles ( |5ikivie & 
Ipser l$9lj). 

The proposals made in this paper can be tested in a va- 
riety of ways. The physical processes can be studied by ap- 
propriate numerical simulations. Conventional iV-body sim- 
ulations of galaxy formation offer limi ted insight into phase 
mixi ng, because of numerical noise ( Hcrna^ms^&^JJtuncj 



old and edge-on they can be used to constrain the disc thick- 
ening and hence angular-momentum transfer rates, (v) Halo 
structure may be directly detectable from small-scale fluc- 
tuations in the surface-brightness distribution of elliptical 
galaxies; these are distinct from surface-brightness fluctu- 
ations arising from bright stars and globular clusters (e.g. 
Jensen, Tonry & Luppino 1998) because they are correlated 
in nearby pixels. Halo structure may also account for dis- 
crepancies in the radio flux ratio s predicted by smooth m ass 
models of gravitational lensing (Mao & Schneider 1996). 



199C 



Hjcrnquist fc Ostriker 199S| ; |5teinmctz fc White 1997TJ 



However, specially designed simulations are more powerful. 
A reasonable approximation is to turn off the gravitational 
forces among the disc particles and among the streamer par- 
ticles, to suppress large-scale instabilities in the disc and 
numerical relaxation within the halo. The only remaining 
forces would be between disc and streamer particles, and 
between both sets of particles and a fixed galactic poten- 
tial. A more complete treatment would include the mutual 
gravitational interactions between streamers. As an illustra- 
tion of the importance of such interactions, consider a single 
streamer consisting of an axisymmetric ring of material in 
an inclined orbit. Eventually this streamer would be dragged 
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